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An exhaustive numerical investigation of the growth of magnetic films in 
confined (d + l)-dimensional stripped geometries {d = 1,2) is carried out by 
means of extensive Monte Carlo simulations. Films in contact with a thermal 
5^ I bath at temperature T, are grown by adding spins having two possible orien- 

tations and considering ferromagnetic (nearest-neighbor) interactions. At low 
temperatures, thin films of thickness L are constituted by a sequence of well- 
ordered domains of average length ^ L. These domains have opposite 
magnetization. So, the films exhibit "spontaneous magnetization reversal" 
during the growth process. Such reversal occurs within a short characteristic 
length Iji, such that Id ^ l-R ^ L. Furthermore, it is found that for d = 1 
the system is non-critical, while a continuous order-disorder phase transition 
at finite temperature takes place in the d = 2 case. Using standard finite-size 
scaling procedures, the critical temperature and some relevant critical expo- 
nents are determined. Finally, the growth of magnetic films in (2 + 1) dimen- 
sions with competing short-range magnetic fields acting along the confinement 
walls is studied. Due to the antisymmetric condition considered, an interface 
between domains with spins having opposite orientation develops along the 
growing direction. Such an interface undergoes a localization-delocalization 
transition that is the precursor of a wetting transition in the thermodynamic 
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limit. Furthermore, the growing interface also undergoes morphological transi- 
tions in the growth mode. A comparison between the well-studied equilibrium 
Ising model and the studied irreversible magnetic growth model is performed 
throughout. Although valuable analogies are encountered, it is found that the 
nonequilibrium nature of the latter introduces new and rich physical features 
of interest. 



1 Introduction 

The preparation and characterization of magnetic nanowires and films is of great 
interest for the development of advanced microelectronic devices. Therefore, the 
study of the behavior of magnetic materials in confined geometries, e.g. thin films, 
has attracted both experimental [jl|, |^, |^, |, |, |], ^, || and theoretical p, [1^, |Tl 



12| , [131 attention. From the theoretical point of view, most of the work has been 



devoted to the study of equilibrium properties of thin magnetic films [^, 10, 11, 12, 



T3| , p!4| and magnetic materials, see e.g. |115| , p!6| , \LTl [181 . In contrast, the aim of 



this work is to study the properties of thin magnetic film growth under far-from- 
equilibrium conditions, by means of extensive Monte Carlo simulations. Within this 
context, this work is related to many recent investigations concerned with irreversible 
growth processes. Indeed, the study of growth systems under far-from-equilibrium 
conditions is a subject that has drawn great attention during the last decades. 
Nowadays, this interdisciplinary field has experienced a rapid progress due to its 
interest in many subfields of physics, chemistry, and even biology, as well as by 
its relevance in numerous technological applications such as the development of 
nanoscale devices, polymer science, crystal and polycrystalline growth, gelation, 
fracture propagation, epidemic spreading, colloids, etc. |^ |21|, |22|, |23 |. 



For the purpose of studying the properties of thin magnetic film growth under 
far-from-equilibrium conditions, a variant of the irreversible Eden growth model pi 



in which particles are replaced by spins that can adopt two different orientations, 
is investigated. Our study is performed in confined (stripped) geometries, which 
resemble recent experiments where the growth of quasi-one-dimensional strips of Fe 
on a Cu(lll) vicinal surface and Fe on a W(llO) stepped substratum [^] have 
been performed. Also, in a related context, the study of the growth of metallic 
multilayers has shown a rich variety of new physical features. Particularly, the 
growth of magnetic layers of Ni and Co separated by a Cu spacer layer has recently 



been studied P5 |. 



The growth of magnetic films with competing short-range magnetic fields, which 
account for the interaction of the growing films with the substrate, is also investi- 
gated in the present work. The competing situation considered leads to rich and 
complex physical phenomena that exhibit a delicate and subtle interplay between 
finite-size effects, wetting, and interface growth mechanisms. Besides the technologi- 
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cal and scientific interest that may arise from tlie nonequilibrium nature of the model 
investigated, as already pointed out, this kind of nonequilibrium wetting phenom- 
ena also appears closely related to very interesting equilibrium wetting transitions, 
which have attracted so far considerable experimental and theoretical attention. For 
instance, surface enrichment and wetting layers have been observed experimentally 
in a great variety of systems, such as e.g. polymer mixtures p6|, |2^, adsorption 
of simple gases on alkali metal surfaces 31], hydrocarbons on mica H, etc. 

From the theoretical point of view, the study of wetting transitions at interfaces has 
been carried out by means of different approaches, such as the mean field Ginzburg- 
Landau method ||3^, transfer matrix and Pfaffian techniques [p5| , p^j , density 
matrix renormalization group methods |3^ , solving the Cahn-Hilliard equation |38| , 
using Molecular Dynamic simulations PU| , solving self-consistent field equations PD 
and by means of extensive Monte Carlo simulations 



41, 42, 43, 



Finally, it should also be remarked that, although the discussion is presented here 
in terms of a magnetic language, the relevant physical concepts could be extended 
to other systems such as fluids, polymers, and binary mixtures. 

This manuscript is organized as follows: in Section 2 details on the model and 
the simulation method are given. Section 3 is devoted to the order- disorder critical 
behavior of magnetic Eden films. Section 4 deals with the study of interfacial phase 
transitions that arise when competing short-range magnetic fields are considered, 
while the conclusions are finally stated in Section 5. 



2 The model and the simulation method 



In the classical Eden model ||2^ on the square lattice, the growth process starts by 
adding particles to the immediate neighborhood (the perimeter) of a seed particle. 
Subsequently, particles are stuck at random to perimeter sites. This growth process 



leads to the formation of compact clusters with a self-affine interface ||2y, ^ |23 
The magnetic Eden model (MEM) ^ considers an additional degree of freedom 
due to the spin of the growing particles. In the present work the MEM is studied in 



{d + 1)— dimensional rectangular geometries for c? = 1, 2, as described in 



For the case d = 1, the MEM is investigated on the square lattice using a 
rectangular geometry L x M (with M ^ L). The location of each site on the lattice 
is specified through its rectangular coordinates (i, j), {X<i<L,l<j< M). The 
starting seed for the growing cluster is a column of parallel-oriented spins placed 
at j = 1 and film growth takes place along the positive longitudinal direction (i.e. 
i > 2). Periodic boundary conditions are adopted along the transverse direction. 
Then, assuming that each spin Sij can be either up or down (i.e. Sij = ±1), clusters 
are grown by selectively adding spins to perimeter sites, which are defined as the 
nearest-neighbor (NN) empty sites of the already occupied ones. Considering a 
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ferromagnetic interaction of strength J > between NN spins, the energy E of a 
given configuration of spins is given by 



I .1 



2 ^ -i-j-i] 

where means a summation taken over all occupied NN sites. 

Analogously, the MEM in (2 + 1) dimensions is studied using a L x L x M 
rectangular geometry (M ^ L). Each site on the lattice is now identified through 
the rectangular coordinates {i,j,k), (1 <i,j<L, l<k< M), and the starting 
seed for the growing film is taken to be a plane of LxL parallel-oriented spins placed 
at /c = 1. The energy of the spin configuration is now calculated by extending the 
summation in Eq.(l) to the three coordinates {i,j,k). For further details on the 
MEM defined in retangular geometries see also [ ^ . 

The last part of this paper (see Section 4) is devoted to the study of the (2 + 1)- 
dimensional MEM with competing short-range magnetic fields applied along one of 
the transverse directions. In this case, the periodic boundary conditions along the 
j— direction are changed to open ones, and competing surface magnetic fields H > 
{H' = —H) acting on the sites placed at j = 1 (j = L) are considered Then, 
the energy of a given configuration of spins in this case will be given by 
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( \ ( 

SijkSi'j'k' ^ ^ \ Silk - Y SiLk ) , (2) 

\{ijk,i'j'k') J (ifc.Si,) 



where {ijk,i j k') means that the summation in the first term is taken over all 
occupied NN sites, while {ik, Si), {ik, T^l) denote summations carried over occupied 
sites on the surfaces Si, (defined as the j = I and j = L planes, respectively). 
Throughout this work we set the Boltzmann constant equal to unity {ks = 1) and 
consider the absolute temperature, energy, and magnetic fields measured in units of 
the coupling constant J. 

The growth process of a MEM film consists in adding further spins to the growing 
film taking into account the corresponding interaction energies. A spin is added to 
the film with a probability proportional to the Boltzmann factor exp{—AE/T), 
where AE is the total energy change involved. At each step, the probabilities of 
adding up and down spins to a given site have to be evaluated for all perimeter sites. 
After proper normalization of the probabilities, the growing site and the orientation 
of the spin are determined through standard Monte Carlo techniques. Although 
both the interaction energy (as given either by Eq.(l) or Eq.(2)) and the Boltzmann 
probability distribution considered for the MEM are similar to those used for the 
Ising model [^], it must be stressed that these two models operate under extremely 
different conditions, namely the MEM describes the irreversible growth of a magnetic 
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material and the Ising model is suitable for the study of a magnetic system under 
equilibrium conditions. In the MEM, the position and orientation of all deposited 
spins remain fixed. During the growth process, the system develops a rough growing 
interface and evolves mainly along the longitudinal direction. Some lattice sites can 
remain empty even well within the system's bulk, but, since at each growth step 
all perimeter sites arc candidates for becoming occupied, these holes arc gradually 
filled. Hence, far behind the active growing interface, the system is compact and 
frozen. When the growing interface is close to reaching the limit of the sample, 
the relevant properties of the irreversibly frozen cluster's bulk (in the region where 
the growing process has definitively stopped) are computed, the useless frozen bulk 
is thereafter erased, and finally the growing interface is shifted toward the lowest 
possible longitudinal coordinate. Hence, repeatedly applying this procedure, the 
growth process is not hmited by the lattice length M. In the present work clusters 
having up to 10^ spins have typically been grown. 

3 Order-disorder critical behavior of magnetic 
Eden films 

Magnetic Eden films grown on a stripped geometry of finite linear dimension L at low 
temperatures show an intriguing behavior that we call spontaneous magnetization 
reversal. In fact, we have observed that long clusters are constituted by a sequence of 
well ordered magnetic domains of average length Ijj^ L. Let Ir be the characteristic 
length for the occurrence of the spontaneous magnetization reversal. Since Ir ~ L, 
we then conclude that the phenomenon has two characteristic length scales, namely 
Ijj and //J, such that Id ^ Ir ^ L. Hence, the spontaneous magnetization reversal 
is essentially due to the small size of the thin film and it becomes irrelevant in 
the thermodynamic limit. Figure 1 shows a snapshot configuration of the (1 + 
1)— dimensional MEM where this phenomenon can be recognized. Here the reversal 
occurring between a domain of spins up (on the left side) and other one constituted 
by spins down (on the right), as well as the interface between both domains, can 
be clearly observed. The magnetization change occurs quite abruptly within the 
characteristic length Ir ~ L. In ordinary thermally driven phase transitions, the 
system changes from a disordered state at high temperatures to a spontaneously 
ordered state at temperatures below some critical value T^, where a second-order 
phase transition takes place. Regarding the Ising model, one has that, in the absence 
of an external magnetic field {H = 0), the low-temperature ordered phase is a 
state with non- vanishing spontaneous magnetization (iM^p). This spontaneous 
symmetry breaking is possible in the thermodynamic limit only. In fact, it is found 
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Figure 1: Spontaneous magnetization reversal observed for L = 32 and T = 0.26 in 
the (1 + 1)— dimensional magnetic film. The snapshot configuration shows the collec- 
tive orientation change: the left (right) domain is constituted by up (down) oriented 
spins. The snapshot corresponds to the bulk of the sample and the growing inter- 
face is not shown. The characteristic length for the occurrence of the magnetization 
reversal, Ir, is of the order of the lattice width, as marked in the figure. 

that the magnetization M of a finite sample formed by N particles, defined by 



N 



M{T,H = 0) = -J2S^iT,H = 0) 



(3) 



can pass with a finite probability from a value near +Msp to another near —Msp, 
as well as in the opposite direction. Consequently, the magnetization of a finite 
system, averaged over a sufficiently large observation time, vanishes irrespective of 
the temperature. That is, the equation M(T, if = 0) ~ holds if the observation 
time {tobs) becomes larger than the ergodic time [terg), which is defined as the time 
needed to observe the system passing from ±Msp to ^Mgp. Since Monte Carlo simu- 
lations are restricted to finite samples, the standard procedure to avoid the problems 
treated in the foregoing discussion is to consider the absolute magnetization as an 
order parameter [Q. Turning back to the MEM, we find that the phenomenon of 
magnetization reversal also causes the magnetization of the whole cluster to vanish 
even for very low (but non-zero) temperatures, provided that the film's total length 
If (which plays the role of tobs) is much larger than Id (which plays the role of terg)- 
Therefore, as in the case of the Ising model [Q, in order to overcome shortcomings 
derived from the finite-size nature of Monte Carlo simulations we have measured the 
mean absolute column magnetization, given by 



\m{j,L,T)\ 



(4) 
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Figure 2: Data corresponding to the (1 + 1)— dimensional MEM: plots of the proba- 
bility distribution of the mean column magnetization PL(m) versus m for the fixed 
lattice width L = 128 and different temperatures, as indicated in the figure. The 
sharp peaks at m = ±1 for T = 0.45 have been truncated, in order to allow a de- 
tailed observation of the plots corresponding to higher temperatures. This behavior 
resembles that of the one-dimensional Ising model. 



It is found that |m(j;L, T)| exhibits a transient growing period with a charac- 
teristic length of order L, followed by the attainment of a stationary regime, which 
is independent of the orientation of the seed. 

The mean column magnetization given by Eq.(4) is a fluctuating quantity that 
can assume L + 1 values. Then, for given values of both L and T, the probability 
distribution of the mean column magnetization {PL{m)) can straightforwardly be 
evaluated, since it represents the normalized histogram of m taken over a sufficiently 
large number of columns in the stationary region [^, ^ . In the thermodynamic 
limit [L — »• oo) the probability distribution {Poo{m)) of the order parameter of an 
equilibrium system at criticality is universal (up to re-scaling of the order parameter) 
and thus it contains very useful and interesting information on the universality class 
of the system For example, Pl{iti) contains information about all 

momenta of the order parameter m, including universal ratios such as the Binder 
cumulant ^2 . 



Figure 2 shows the thermal dependence of Piini) versus m, as obtained for the 
(1 + 1)— dimensional MEM. We can observe that at high temperatures PL{m) is 
a Gaussian centered at m = 0. As the temperature gets lowered, the distribution 
broadens and develops two peaks at m = ±1. Further decreasing the temperature 
causes these peaks to become dominant while the distribution turns distinctly non- 
Gaussian, exhibiting a minimum just at m = 0. It should be pointed out that the 
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emergence of the maxima at m = ±1 is quite abrupt. This behavior reminds us the 
order parameter probabihty distribution characteristic of the one-dimensional Ising 
model. In fact, for the well studied d— dimensional Ising model we know 

that for T > T^, Pz^(M) is a Gaussian centered at M = 0, given by 

Pz.(M) (X exp (^^^j , (5) 

where the susceptibility x is related to order parameter fluctuations by 

Decreasing the temperature, the order parameter probability distribution broadens, 
it becomes non-Gaussian, and near Tc it splits into two peaks that get the more sep- 
arated the lower the temperature. For T < Tc and linear dimensions L much larger 
than the correlation length ^ of order parameter fluctuations, one may approximate 
Pl{M) near the peaks by a double-Gaussian distribution, i.e. 

Pl{M) cx exp y ^ j + exp ^ ^ j , (7) 

where Mgp is the spontaneous magnetization, while the susceptibility x is now given 
by 

X= y((M2)-(|M|)^) . (8) 

From Eq.(5) it turns out that the Gaussian squared width associated with high 
temperature distributions is very close to the 2nd moment of the order parameter, 
i.e. 



a 



2 ^ / » ,r2 



(M^) . (9) 



It should be noticed that this equation is a straightforward consequence of the 
Gaussian shape of the order parameter probability distribution and, thus, it holds 
for the MEM as well. From the well known one-dimensional exact solution for a 
chain of L spins one can establish the relationship 



X=^exp(2/T) ; (10) 

then, Eqs.(6) and (10) lead us to 

(M2) = iexp(2/T) (11) 

(where it has been taken into account that (M) = due to finite-size effects, irrespec- 
tive of temperature). From Eqs.(9) and (11) we can see that the high-temperature 
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Figure 3: Data corresponding to the (2 + 1)— dimensional MEM: plots of the prob- 
ability distribution PL{m) versus m for the fixed lattice size L — 16 and different 
temperatures, as indicated in the figure. The occurrence of two maxima located at 
m = ±Msp (for a given value of Mgp such that < Mgp < 1) is the hallmark of a 
thermal continuous phase transition that takes place at a finite critical temperature. 

Gaussian probability distribution broadens exponentially as T gets lowered, until it 
develops delta-like peaks at M = ±1 as a consequence of a boundary effect on the 
widely extended distribution. It should be noted that for d > 2 this phenomenon 
is prevented by the finite critical temperature which splits the Gaussian, as implied 
byEq.(7). 

Figure 3 shows the thermal evolution of the probability distribution as obtained 
for the (2 -|- 1)— dimensional MEM. Notice that now m is defined by an average 
over transverse planes constituted by L x L spins, analogously to Eq.(4). Hence, m 
takes now + 1 possible values. For high temperatures, the probability distribution 
corresponds to a Gaussian centered at m = 0. At lower temperatures we observe 
the onset of two maxima located at m = ±Msp (0 < Msp < 1), which become 
sharper and approach m = ±1 as T is gradually decreased. These low-temperature 
probability distributions clearly reflect the occurrence of the magnetization reversal 
effect already discussed for the case of (1 + 1)— dimensional magnetic fllms. 

Figure 4 shows the location of the maximum of the probability distribution as a 
function of temperature for both (d + l)-dimensional MEM models (with d = 1,2), 
where only maxima located at m > are considered, since the distributions are 
symmetric around m = 0. After inspection of figure 4, it becomes apparent the 
different qualitative behavior of both systems. Indeed, while for the d = 2 case 
we observe a smooth transition from the rrimax = value characteristic of high 
temperatures to nonzero rrimax values that correspond to lower temperatures, the 
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Figure 4: Plots showing the location of the maximum of the probability distribution 
as a function of temperature for both {d + 1)— dimensional MEM models {d = 1,2). 
The lines are guides to the eye. The smooth transition for d=2 is the signature of 
a thermal continuous phase transition occurring at a finite critical temperature. 



curve obtained for d = 1 shows instead a Heaviside-like jump. In contrast to the 
(1 + 1)— dimensional case, the behavior exhibited by the (2 + l)-dimensional MEM 
(e.g. as displayed by figures 3 and 4) is the signature of a thermal continuous phase 
transition that takes place at a finite critical temperature. 

From the finite-size scaling theory, developed for the treatment of finite-size ef- 



fects at criticality and under equilibrium conditions |5g, it is well known that 
if a thermally driven phase transition occurs at a temperature Tc > in the ther- 
modynamic limit, then in a confined geometry of linear dimension L this transition 
becomes smeared out over the temperature region AT(L) around a shifted effec- 
tive transition temperature Tc{L), which are related to L through phenomenological 
exponents. Indeed, it is found that 

AT(L) oc L-^ (12) 

and 

|r,(L) - r,| cx L-^ , (13) 

where the rounding and shift exponents are given hj 6 = X = u^^, respectively, and 
where u is the exponent that characterizes the divergence of the correlation length 
at criticality. 

Furthermore, from well-established finite-size scaling relations, the following 
Ansatze hold just at criticality: 

(|m(L,T = T,)|)ocL-^/^ (14) 

and 

XmaAL) ^ L-^^" , (15) 
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Figure 5: Plots of the effective finite-size critical temperatures Tcn{L) versus (for 
n— 1,2) corresponding to the (2 + 1)— dimensional magnetic film. Tci(L) is defined 
as the value that corresponds to (|m|) = 0.5, while Tc2{L) is the temperature that 
corresponds to the maximum of the susceptibility. The solid lines show the linear 
extrapolations that meet at the critical point given by Tc — 0.69 ± 0.01. 



where (3 and 7 are the order parameter and the susceptibility critical exponents, 
respectively. Note that Xm,ax{L), as given by Eq.(15), refers to the maximum of 
x{L,T) as a function of T for a fixed lattice size L. 

In order to describe quantitatively the critical behavior of the MEM in (2 + 1)- 
dimensions, we may test the validity of the scaling relations given in Eqs.(12)-(15). 
As in the case of equilibrium systems, in the present case various "effective" L- 
dependent critical temperatures can also be defined. In particular, we will define 
Tci{L) as the value that corresponds to (|m|) = 0.5 for fixed L, and Tc2{L) as the 
one corresponding to the maximum of the susceptibility for a given L, assuming that 
the susceptibility is related to order parameter fluctuations in the same manner as 
for equilibrium systems (as given by Eqs.(6) and (8)). Then, Tc can be obtained 
from plots of Tc„(L) versus L^^ (for n = 1, 2), as is shown in figure 5. Following this 
procedure we find that both Tci{L) and Tc2{L) extrapolate (approximately) to the 
same value, allowing us to evaluate the critical temperature Tc — 0.69 ± 0.01 in the 
thermodynamic limit. 

After determining Tc, the correlation length exponent v can be evaluated by 
means of Eq.(13), making the replacement \ = Indeed, taking Tc at the mean, 
maximum and minimum values allowed by the error bars, we obtain six log-log plots 
of |Tc„(L) — Tc| versus L for n = 1, 2. The slope of each of these plots, not shown 
here for the sake of space, yields a value for v. The obtained values are: 

i/ = 1.08 (Tc = 0.68), v^im (Tc = 0.69), 
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Figure 6: Log-log plot of Xmax versus L. The linear fit (solid line) yields 
7/1/ = 2.02 ± 0.04. 

i/ = 0.88 (Te^O.TO) for n = 1, (16) 

and 

i/ = 1.20 (7; = 0.68), i/ = 1.08 (7; = 0.69), 

i/ = 0.95 (7; = 0.70) for n = 2. (17) 

Thus our estimate is given by = 1.04 ±0.16, where the error bars reflect the error 
derived from the evaluation of Tc, as well as the statistical error. 

Studying the susceptibility x as a function of the temperature for several different 
lattice sizes, it is found that x exhibits a peak, which becomes sharper and shifts 
toward lower temperatures as L is increased. Hence, Eq.(15) can be used to evaluate 
7/1/ from the slope of a log-log plot of Xmax versus L, as figure 6 shows. The linear 
fit yields ^/v — 2.02 ± 0.04. Using this value and the value formerly obtained for v 
we thus determine 7 = 2.10±0.36. 

Figure 7 shows log-log plots of {\m\){T = T^) versus L for the mean, maximum 
and minimum allowed values of Tc. Considering only the larger lattices, the linear 
fits to the data according to Eq.(14) yield the following estimates: P/u — 0.11, 
P/u = 0.16 and P/iy = 0.19. We then assume the value P/u = 0.15 ± 0.04, where 
the error bars refiect the error derived from the evaluation of Tc, as well as the 
statistical error. From this value and the value formerly obtained for u we thus 
determine /3 = 0.16 ± 0.05. 

In this manner, we conclude that magnetic Eden films grown on stripped {d+1)- 
dimensional geometries are non-critical for d = 1, in the sense that the ordered phase 
is trivially found only at T = 0. However, MEM films exhibit a continuous order- 
disorder phase transition at the critical temperature Tc — 0.69 ± 0.01 for d — 2. 



12 




Figure 7: Log-log plots of average magnetization measured at the critical point 
(m(T = Tc)) versus as obtained for the mean, maximum and minimum allowed 
values of Tg. The hnear fits (sohd lines) yield an estimate 13 /v — 0.16 ± 0.05. 

Furthermore, the critical exponents of the MEM in (2 + 1) dimensions, as obtained 
using a finite-size scaling analysis, are: u = 1.04 ± 0.16, 7 = 2.10 ± 0.36, and 
(3 = 0.16 ±0.05. 

4 Study of the (2 + l)-dimensional MEM with 
competing short-range magnetic fields: Wet- 
ting and morphological phase transitions 

In this section, we will study the interfacial phase transitions that arise in the (2-1-1)- 
dimensional MEM, when competing short-range magnetic fields applied along one of 
the transverse directions are considered. As described in Section II, we will assume 
competing surface magnetic fields H > {H' — —H) acting on the surfaces j — 1 
and j = L. Hence, the energy associated to a given spin configuration acquires an 
additional term due to the interaction of the surface spins with the applied magnetic 
fields (see Eq.(2)). 

As shown in figure 8, magnetic Eden films that grow in a confined geometry with 
competing surface fields exhibit a very rich phase diagram, which is composed of 
eight regions. These regions are delimited by several distinct, well-defined transition 
curves. As will be shown below, the bulk order-disorder (finite-size) critical point 
Tc(L), the Ising-like quasi- wetting transition curve Tyj{L,H), and two morpholog- 
ical transitions associated to the curvature of the growing interface (namely, from 
convex to non-defined to concave), can be quantitatively located. Moreover, in or- 
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Figure 8: H — T phase diagram corresponding to a lattice of size L = 12. The verti- 
cal straight line at Tc{L) = 0.84 corresponds to the L— dependent critical tempera- 
ture, which separates the low-temperature ordered phase from the high-temperature 
disordered phase. Open (filled) circles refer to the transition between non-defined 
and concave (convex) growth regimes, and squares stand for the Ising-like local- 
ization-delocalization transition curve. Eight different regions are distinguished, 
as indicated in the figure. Also indicated are seven representative points that are 
discussed in the text. The inset shows the phase diagram corresponding to the 
thermodynamic limit composed of six different regions. 

der to gain some insight into the physics involved in this complex phase diagram, 
some typical snapshot configurations characteristic of the various different growth 
regimes are obtained (see figure 9) and discussed. Finally, the phase diagram in the 
thermodynamic limit will be drawn (see inset of figure 8) by extrapolating finite-size 
results (see figure 10). 

The (L— dependent) bulk order-disorder critical temperature can be identified 
with the peak of the susceptibility at zero surface field. For L = 12, the critical point 
so defined is Tc{L = 12) = 0.84, and is shown in figure 8 by a vertical straight line. 
So, the left (right) hand side part of the phase diagram corresponds to the ordered 
(disordered) growth regime that involves Regions /, //, ///, IV, and A (Regions 
V,VI, and B). 

Using a standard procedure §], the localization-delocalization transition curve 
(on the H — T plane) corresponding to the up-down interface running along the walls 
can be computed, considering that a point with coordinates {H^,T^) on this curve 
maximizes x{H,T). So, the size-dependent localization-delocalization transition 
curve is obtained, as shown in figure 8 (open squares). As in the case of the Ising 
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model, this quasi-wetting transition refers to a transition between a nonwet state 
that corresponds to a locahzed interface bound to one of the confinement walls, and 
a wet state associated with a delocalized domain interface centered between roughly 
equal domains of up and down spins. The localization-delocalization transition in a 
confined system is indeed the precursor of the true wetting transition, which occurs 
in the thermodynamic limit ^ In fact, there is observed a finite jump in the 
wetting layer thickness that takes place as a result of the finite size of the system. 
As the lattice size is increased, the magnitude of the jump grows and diverges in the 
L ^ oo limit, as expected for a continuous wetting transition. 

Since the MEM is a nonequilibrium kinetic growth model, it also allows the 
identification of another kind of phase transition, namely a morphological transi- 
tion associated with the curvature of the growing interface of the system. To avoid 
confusion, we would like to remark that the term interface is used here for the trans- 
verse interface between occupied and empty lattice sites, while it was used above 
for the longitudinal interface between up and down spin domains. To explore this 
phenomenon quantitatively, the behavior of the contact angles between the growth 
interface and the confinement walls (as functions of temperature and magnetic field) 
have to be investigated. Clearly, two different contact angles should be defined in 
order to locate this transition, namely 6d for the angle corresponding to the domi- 
nant spin cluster, and 9nd for the one that corresponds to the non-dominant spin 
cluster. Both contact angles can straightforwardly be determined by measuring the 
location of the growth interface averaged over a sufficiently long growing time. In 
this way, three different growth regimes can be distinguished: (i) the concave growth 
regime that occurs when the system partially wets the walls on both sides (i.e. for 
Gdi^nd < f), (ii) the convex growth regime that occurs for 6e),6nd > f , and (Hi) 
the regime of non-defined curvature that occurs otherwise. The corresponding tran- 
sition curves obtained for confined magnetic Eden films are shown on the H — T 
phase diagram of figure 8, where open (filled) circles refer to the transition between 
non-defined and concave (convex) growth regimes. 

As anticipated above, we will now introduce and discuss some characteristic 
snapshot pictures, in order to provide qualitative explanations that account for the 
different growth regimes observed. Let us begin with Region I (see figure 8), that 
corresponds to the Ising-like nonwet state and the convex growth regime. In this 
region, the temperature is low and the system grows in an ordered state, i.e. the 
dominant spin domain prevails and the deposited particles tend to have their spins 
all pointing in the same direction. Small clusters with the opposite orientation may 
appear preferably on the surface where the non-dominant orientation field is applied. 
These "drops" might grow and drive a magnetization reversal, thus changing the 
sign of the dominant domain. Indeed, the formation of sequences of well-ordered 
domains are characteristic of the ordered phase of confined (finite-size) spin systems 
such as the Ising magnet [Q. Due to the open boundary conditions, perimeter sites 
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Figure 9: Snapshot pictures showing a longitudinal slice given by a fixed value of 
the transverse coordinate i. Grey (black) points correspond to spins up (down). 
The surface field on the upper (lower) confinement wall is positive (negative). The 
snapshots correspond to a lattice size L = 32 and several different values of temper- 
ature and surface fields: {&)H = 0.05, T = 0.6; {h)H = 0.5, T = 0.55; {c)H = 1.4, 
T = 0.6; {d)H = 0.1, T = 1.0; {e)H = 1.6, T = 1.4; and (f) H = 0.20, T = 0.82. 



at the confinement walls experience a missing neighbor effect, that is, the number 
of NN sites is lower than for the case of perimeter sites on the bulk. Since H is too 
weak to compensate this effect, the system grows preferentially along the center of 
the sample as compared to the walls, and the resulting growing interface exhibits 
a convex shape. So, Region / corresponds to the Ising-like nonwet state and the 
convex growth regime. A typical snapshot configuration characteristic of Region / 
is shown in figure 9(a). 

Let us now consider an increase in the surface magnetic fields, such as the sys- 
tem may be driven into Region II (see figure 8). Since the temperature is kept 
low, the system is still in its ordered phase and neighboring spins grow preferably 
parallel-oriented. The surface fields in this region are stronger and thus capable of 
compensating the missing NN sites on the surfaces. But, since the fields on both 
surfaces have opposite signs, it is found that, on the one hand, the field that has 
the same orientation as the dominant spin cluster favors the growth of surface spins, 
while on the other hand, the sites on the surface with opposite field have a lower 
probability to be chosen during the growth process. Hence, the contact angle cor- 
responding to the dominant spin cluster is then ^z? < f , while the non-dominant is 
Ond > f • Thus, on the disfavored side the growing interface becomes pinned and 
the curvature of the growing interface is not defined. Figure 9(b) shows a typical 
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snapshot corresponding to Region //. 

Keeping H fixed within Region // but increasing the temperature, thermal noise 
will enable the formation of drops on the disfavored side that eventually may nucle- 
ate into larger clusters as the temperature is increased even further. This process 
may lead to the emergence of an up-down interface, separating oppositely oriented 
domains, running in the longitudinal direction (i.e. parallel to the walls). Since sites 
along the up-down interface are surrounded by oppositely oriented NN spins, they 
have a low growing probability. So, in this case the system grows preferably along 
the confinement walls and the growing interface is concave (figure 9(c)). Then, as 
the temperature is increased, the system crosses to Region A (see figure 8) and the 
onset of two competitive growth regimes is observed, namely: (i) one exhibiting a 
non-defined growing curvature that appears when a dominant spin orientation is 
present, as in the case shown in figure 9(b); (ii) another that appears when an up- 
down interface is established and the system has a concave growth interface, as is 
shown in figure 9(c). 

On further increasing the temperature and for large enough fields, the formation 
of a stable longitudinal up-down interface that pushes back the growing interface 
is observed. So, the system adopts the concave growth regime (see figure 9(c) 
corresponding to Region IV in figure 8). Increasing the temperature beyond Tc{L), a 
transition from a low-temperature ordered state (Region IV) to a high-tcmpcrature 
disordered state (Region VI, see figure 9(e)), both within the concave growth regime, 
is observed. Analogously, for small enough fields, a temperature increase drives the 
system from the ordered convex growth regime (Region I) to the disordered convex 
growth regime (Region V, see figure 9(d)). As shown in figure 8, there is also an 
intermediate fiuctuating state (Region B) between Regions V and VI, characterized 
by the competition between the disordered convex growth regime and the disordered 
concave one. 

Finally, a quite unstable and small region (Region III in figure 8) that exhibits 
the interplay among the growth regimes of the contiguous regions, can also be iden- 
tified. Since the width of Region III is of the order of the rounding observed in 
Tc{L), large fiuctuations between ordered and disordered states are observed, as well 
as from growth regimes of non-defined curvature to convex ones. However, figure 
9(f) shows a snapshot configuration that is the fingerprint of Region III, that may 
prevail in the thermodynamic limit, namely a well defined spin up-down interface 
with an almost flat growing interface. 

Let us now extrapolate our results to show that the rich variety of phenomena 
found in a confined geometry is still present in the thermodynamic limit (L — > oo), 
leading to the phase diagram shown in the inset of figure 8. As clearly seen by 
comparison with the finite-size results, the crossover Regions A and B collapse in 
this limit, so only the six regions that correspond to well identified growth regimes 
(as illustrated by the snapshot configurations of figure 9) appear to remain. 
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Figure 10: Plots of T versus for 12 < L < 48, corresponding to the points 
Pi, P*, P2, and P2*) all of them with H — 0.6. The fits to the data (solid hues) show 
that, within error bars, Pj — > P* (i = 1, 2) for L — > 00. 

In order to illustrate the extrapolation procedure, the following seven represen- 
tative points of the finite-size phase diagram are discussed in detail: (i) the points 
labeled Pi, Pj*, P2, and P2*, that correspond to the intersections of the H = 0.6 
line with the various transition curves shown in figure 8, and (ii) the points labeled 
P3, P|, and P4, that refer to the intersection point between Regions /, //, ///, and 
A, the minimum of the limiting curve between Regions IV-VI and A-B, and the 
zero-field transition point, respectively. 

Figure 10 shows plots of T versus for 12 < L < 48 corresponding to the 
points Pi,P*,P2, and Pg*. Also shown in the figure are the fits to the data ex- 
trapolated to — 0. The results from the extrapolations are: Ti = 0.67 ± 0.01, 
T* = 0.66 ± 0.01, and T2 = 1.30 ± 0.02, T* = 1.29 ± 0.01, pointing out that, within 
error bars, Pj P* {i = 1, 2) in the L 00 limit. Using the same procedure, the ex- 
trapolations of P3 and P3* (not shown here) give: H3 = 0.30±0.01, = 0.31±0.02, 
and T3 = 0.69 ± 0.01, T* = 0.71 ± 0.03. So, one has P3 P3* for L ^ 00 within 
error bars. Finally, the extrapolation of P4 is T4 — TC — 0.69 ± 0.01. 

Using the above-mentioned extrapolation procedure, the phase diagram in the 
thermodynamic limit can be drawn, as shown in the inset of figure 8. By comparison 
with the finite-size phase diagram of figure 8, one can note that, as anticipated, the 
crossover Regions A and B appear in the phase diagram just as a consequence of 
the finite-size nature of confined geometries, since they collapse in the L — > 00 limit. 
Moreover, we conjecture that Region /// may remain in the thermodynamic limit. 
Although this (very tiny!) region corresponds to a physically well-characterized 
growth regime, since one expects that the system in this region may grow in an 
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ordered phase with a delocahzed up-down domain interface and a convex growing 
interface, statistical errors due to large fluctuations close to criticality hinder a more 
accurate location of this region. The unambiguous clarification of our conjecture 
remains as an open question that will require a huge computational effort. 

Besides an Ising-like continuous wetting transition, coupled morphological tran- 
sitions in the growing interface, which arise from the MEM's kinetic growth process, 
have also been identified. Comparing the equilibrium wetting phase diagram of the 
Ising model |33|, ^ and that of the MEM, it follows that the nonequilibrium 



nature of the latter introduces new and rich physical features of interest: the non- 
wet (wet) Ising phase splits out into Regions / and II (Regions /// and IV), both 
within the ordered regime (T < Tc) but showing an additional transition in the in- 
terface growth mode. Also, the disordered state of the Ising system (T > Tc) splits 
out into Regions V and VI exhibiting a transition in the interface growth mode. 



5 Conclusions 

In the present work we have studied the growth of magnetic Eden films with fer- 
romagnetic interactions between nearest-neighbor spins in a (d -|- 1)— dimensional 
rectangular geometry (for d = 1,2), by means of extensive Monte Carlo simula- 
tions. For both dimensions the phenomenon of spontaneous magnetization reversal 
is observed at low temperatures. Indeed, MEM films grown at low temperatures are 
constituted by a sequence of magnetic domains, each of them with a well-defined 
magnetization, such that the magnetization of adjacent domains is antiparallel. Fur- 
thermore, it is found that the (1 -|- 1)— dimensional MEM is non-critical, while the 
(2 -|- 1)— dimensional MEM undergoes a thermally driven second-order phase tran- 
sition at finite temperature, which is evaluated by extrapolating some "effective" 
L-dependent critical temperatures to the thermodynamic limit (L oo). Using a 
finite-size scaling theory, some relevant critical exponents that characterize the be- 
havior of the (2 -|- 1)— dimensional MEM at criticality are determined. The observed 
behavior is reminiscent to that of the equilibrium Ising model, although it should 
be stressed that the MEM is a far-from-equilibrium growing system. 

Finally, the (2 + 1)— dimensional MEM with competing surface magnetic fields, 
which may account for the interaction of the growing magnetic films with the sub- 
strate, is investigated. An Ising-like localization-delocalization wetting transition 
and, on the other hand, a morphological transition associated with the curvature of 
the growing interface, are located. In this way, eight different regions on the H — T 
phase diagram for a finite-size lattice are identified. Moreover, the characteristic be- 
havior of typical growth processes within each region are discussed, and qualitative 
explanations that account for the observed features are provided. Extrapolating the 
results obtained for various lattice sizes, the phase diagram corresponding to the 
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L — > oo limit is also determined. It is composed of six different regions, since two 
crossover regions identified in the finite-size phase diagram appear to collapse in the 
thermodynamic limit. The phase diagram obtained shows new and rich physical 
features of interest, which arise as a consequence of the nonequilibrium nature of 
the model investigated. 

We expect that the present stndy will contribute to the fields of irreversible 
growth processes in confined geometries and nonequilibrium wetting phenomena, 
and we hope that it will stimulate further experimental and theoretical work in 
these topics of widespread technological and scientific interest. 
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